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MULTDsEAM  RADAR  ALTIMETRY:  SPACEBORNE  FEASIBILITY 
AND  AIRBORNE  EXPERIMENTATION 


INTRODUCTION 


The  U.S.  Navy  has  a  high-priority  need  for  an  eddy-resolving  global  ocean  prediction  system  [1] 
to  support  antisubmarine  warfare  (ASW)  and  other  naval  operations.  This  system  [2-4]  should  con¬ 
tain  an  adequately  sampled  satellite  data  input,  global  surface-wind  and  heat-flux  inputs  from  atmos¬ 
pheric  models,  a  large  computer  with  a  sustained  speed  of  gigaflops  and  100  million  words  of 
memory,  and  efficient  and  validated  ocean  models  for  data  assimilation  and  forecasting. 

As  demonstrated  by  GEOS  AT  [5],  one  can  infer  subsurface  oceanographic  parameters  with 
dynamic  heights  measured  by  the  altimeter.  This  potential  to  measure  dynamic  heights  and  an  all- 
weather  capability  make  an  altimeter  the  most  promising  source  of  global  input  data  for  future  eddy¬ 
resolving  global  ocean  prediction  systems  [2].  The  main  drawback  of  present-day  satellite  altimetry  is 
that  only  subsatellite,  or  line,  mapping  is  furnished.  To  provide  adequate  spatial  and  temporal  cover¬ 
age,  a  single  multibeam  altimeter  or  a  multisatellite  system  with  single-beam  altimeters  is  needed. 

From  the  ASW  and  acoustic  modeling  point  of  view  [6],  the  most  tactically  significant  eddies 
are  those  with  diameters  that  exceed  100  km.  The  typical  surface  swirl  velocity  of  these  mesoscale 
features  is  in  the  10  cm/s  to  1  m/s  range.  To  estimate  the  precision  and  beam  configuration  that  a 
multibeam  altimeter  should  have,  the  geostrophic  assumption  can  be  used  to  compute  the  sea  surface 
height  anomaly  Ah  that  would  be  generated  by  the  eddies  in  the  open  ocean.  If  C  is  the  swirl 
velocity,  D  is  the  diameter  of  a  circular  eddy,  g  is  the  acceleration  due  to  gravity,  <t>  is  the  geographic 
latitude,  and  Q  is  the  angular  velocity  of  the  earth,  the  geostrophic  analysis  gives 


Ah  ~ 


0  sin  <(> 
8 


CD. 


Figure  1  plots  this  relation.  It  shows  that  an  altimeter  with  a  height  measurement  precision  of  S 
cm  could  detect  eddies  whose  velocity  and  diameter  combination  fall  into  the  area  to  the  right  of  the 
Ah  —  5  cm  curve  at  the  three  indicated  latitudes.  A  precision  of  s5  cm  is  required  to  detect  the 
most  tactically  significant  eddies. 

Figure  2  shows  the  type  of  coverage  a  three-beam  altimeter  has  when  the  beam  spot  separation 
on  the  sea  surface  is  50  km.  The  orbit  shown  is  for  an  800-km  altitude  GEOSAT-type  mission  with  a 
17-day  exact-repeat  track  pattern.  The  solid  lines  denote  the  coverage  if  only  a  single  nadir  beam  is 
used.  The  dotted  lines  denote  the  off-nadir  beams.  The  background  in  Fig.  2  is  a  phytoplankton  pig¬ 
ment  image  produced  by  the  Coastal  Zone  Color  Scanner  aboard  NASA’s  NIMBUS-7  (shown  here  in 
black).  A  warm  core  eddy  is  the  circular  feature  in  the  lower  left  of  the  figure. 
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Fig.  1  —  An  altimeter  with  a  high  measurement  precision 
(Ah)  >5  cm  could  detect  eddies  whose  velocity  and 
diameter  combinations  fall  into  the  area  to  the  right  of  the 
curve  at  three  representative  latitudes 
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Fig.  2  —  The  orbit  tracks  are  for  an  800-km  altitude  GEOSAT-type  orbit  with  a  17-day  exact-repeat  track  pattern. 
The  solid  lines  denote  the  coverage  if  only  a  single  nadir  beam  is  used.  The  dotted  lines  denote  two  off  nadir 
beams.  The  beam  spot  separation  on  the  equator  is  50  km. 


From  a  measurement  point  of  view,  a  single  altimeter  having  three  beams  and  5-cm  precision  or 
three  satellites  each  with  a  single  5-cm  precision  altimeter  are  possible  configurations.  From  a  multi¬ 
beam  altimeter  system  design  point  of  view,  a  number  of  questions  exist  regarding  the  viability  of  the 
multibeam  approach.  Salient  among  these  are  the  inherent  height  and  spatial  resolution  achievable 
with  realistic  systems  and  the  pointing  angle  accuracy.  This  report  focuses  on  these  and  other 
unresolved  areas. 

We  consider  only  the  problem  of  mesoscale  mapping  with  a  beam  ground  track  separation  of  50 
km.  Figure  3  shows  the  typical  geometry.  The  next  section  of  this  report  first  identifies  potential 
system  configurations  for  accomplishing  off-nadir  altimetry  and  then  presents  a  quantitative  analysis 
of  the  candidate  systems.  The  final  section  discusses  aircraft  experiments  that  are  recommended  to 
provide  experimental  corroboration  between  NASA/Wallops  and  the  Naval  Research  Laboratory. 

SPACEBORNE  FEASIBILITY 

In  conventional  nadir-oriented  radar  altimetry,  the  leading  edge  or  ramp  portion  of  the  average 
backscattered  return  is  proportional  in  its  length  to  the  width  of  the  transmitted  pulse.  The  trailing 
edge  of  the  return  is  determined  primarily  by  the  antenna  pattern  and.  to  a  lesser  degree,  the  scatter¬ 
ing  cross  section  of  the  surface  [7 1 .  The  ocean  surface  waveheight  is  measured  by  elongating  or 
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Fig.  3  —  Geometry  for  maximum  cross-track  beam 


stretching  the  leading  edge  of  the  average  return.  The  altitude  measurement  is  also  dependent  on  the 
shape  of  the  average  return.  Figure  4(a)  shows  an  average  return  waveform  for  a  conventional  nadir 
altimeter  operating  over  the  ocean;  Fig.  4(b)  shows  a  typical  off-nadir  return. 

When  an  altimeter  beam  moves  off-nadir,  the  influence  of  the  waveheight  and  the  transmitted 
pulse  shape  ou  the  mean  return  waveform  is  diminished.  In  fact,  when  the  off-nadir  angle  exceeds  a 
very  few  beamwidths,  the  shape  of  the  averaged  waveforms  is  almost  exclusively  determined  by  the 
antenna  beamwidth  [8].  In  this  case,  the  average  waveform  shape  is  essentially  independent  of  the 
ocean  waveheight.  Hence,  measuring  waveheight  does  not  appear  to  be  possible  for  off-nadir  radar 
beams.  However,  the  possibility  still  exists  to  measure  the  range  to  the  center  of  the  off-nadir  beam 
and  the  power  backscattered  to  the  radar  from  this  area  on  the  surface.  These  two  measurements  hold 
the  promise  of  both  wide  swath  topographic  profiling  and  wind  speed  measurement. 

The  above  discussion  is  intended  to  define  those  quantities  that  can  be  measured  from  an  off- 
nadir  attitude  and  those  that  are  measurable  but  with  degraded  precision.  For  these  latter  quantities, 
the  feasibility  of  measuring  these  quantities  is  directly  related  to  precision.  Consequently,  the 
remainder  of  this  report  is  devoted  to  estimating  the  precision  with  which  a  given  technique  can 
measure  the  range  from  the  altimeter  to  the  surface  by  using  off-nadir  pointed  beams.  The  question 
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TIME 

Fig.  4a  —  Average  return  waveform  for  a  conventional  nadir  altimeter 
operating  over  the  ocean 


Fig.  4b  —  Typical  off  nadir  return  plot 
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of  the  radar  cross  section  measurement  is  not  directly  addressed  here  because  this  measurement  is 
inherently  much  less  difficult.  That  is,  the  range,  or  altitude  measurement,  is  significantly  more  com¬ 
plicated  than  the  return  power  measurement. 

CANDIDATE  SYSTEMS 

This  section  describes  proposed  means  for  accomplishing  wide  swath,  or  off- nadir,  altimetry. 
Subsequent  sections  consider  the  precision  of  these  techniques. 

Pure  Multiple  Beam 

Brown  [9]  was  the  first  to  propose  multiple-beam  altimetry  to  measure  cross-track  surface  slope. 
This  system  used  a  five-beam  configuration  comprising  a  nadir  directed  beam,  right-  and  left-of-track 
(cross-track)  beams,  and  fore-  and  aft-pointed  (along-track)  beams.  The  two  cross-track  beams  and 
the  nadir  beam  measured  the  cross-track  slope,  and  the  fore-  and  aft-pointed  beams,  in  conjunction 
with  the  nadir  beam,  measured  along  track.  The  along-track  measurement  could  be  compared  with 
the  continuous  recording  by  the  nadir  beam  to  estimate  the  precision  of  the  three-beam  scheme  and  to 
provide  precursor  information. 

Brown  [9]  did  not  estimate  the  precision  of  this  slope-measuring  seneme.  Since  the  off-nadir 
beams  pointing  angles  were  restricted  to  no  more  than  eight  beamwidths,  the  antenna  for  this  system 
was  an  offset-fed  parabola.  For  very  large  antennas,  the  pointing  error  was  recognized  as  a  potential 
problem,  but  using  the  nadir  beam  to  aid  in  the  attitude  determination  was  considered  a  practical  solu¬ 
tion.  This  system  is  called  a  pure  multiple-beam  system  because  it  uses  one  independent  beam  per 
measurement  cell.  Of  all  the  systems  discussed,  this  is  the  most  straightforward  in  concept  and  easi¬ 
est  to  implement. 

Interferometric  Approach 

A  drawback  of  the  proposed  pure  multiple-beam  system  is  that  it  relies  on  very  narrow  beams  to 
achieve  the  off-nadir  topographic  measurements.  Bush  et  al.  [10]  proposed  the  use  of  a  two-element 
interferometer  to  develop  very  narrow  beams  along  with  phasing  the  elements  to  shift  the  beams  off 
nadir.  In  this  approach,  the  width  of  the  lobes  ic  the  interferometer  radiation  pattern  is  determined  by 
the  spacing  between  the  two  antenna  elements  of  the  array.  The  pattern  of  each  element  attenuates 
the  array  lobes  adjacent  to  the  desired  lobe.  McGoogan  and  Walsh  [11]  considered  the  two-element 
interferometer  to  be  more  suitable  for  space  application  than  the  pure  multiple-beam  system.  In  1985, 
the  Applied  Physics  Laboratory  of  Johns  Hopkins  University  [12]  developed  a  conceptual  design  of  an 
interferometric  system  for  space  borne  applications.  This  system  used  two  1.5-m  parabolic  reflectors 
displaced  by  5  m  left  of  tracK  and  5  m  right  of  track  for  a  total  interferometer  baseline  distance  of  10 
m.  With  a  50-km  one-sided  beam  squint,  the  precision  of  the  resulting  range  measurement  was 
predicted  to  be  5  cm  with  a  5  kHz  pulse  repetition  frequency  (PRF). 

Synthetic  Aperture  Radar  Approach 

Two  categories  of  synthetic  aperture  radar  (SAR)  are  available  for  use  with  global  oceano¬ 
graphic  altimetry.  First,  in  its  simplest  implementation,  a  real  aperture  can  provide  cross-track  reso¬ 
lution  and  a  synthetic  aperture  (or  Doppler  beam  sharpening  [13])  can  provide  along-track  resolution 
[14].  With  this  approach,  a  narrow  rectangular,  or  fan-beam,  antenna  would  be  used.  In  fact,  a 
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quasi-optics  antenna  would  be  required,  since  a  traveling-wave  antenna  would  radiate  incrementally  as 
the  short  pulse  ( ~2  ns)  propagated  along  the  aperture. 

In  the  second  SAR  category,  the  Naval  Research  Laboratory  (NRL)  and  others  [14,15]  have 
investigated  the  use  of  SAR  systems  within  an  interferometer  configuration.  In  principle,  this  tech¬ 
nique  can  provide  contiguous  cross-track  altimetry  information.  However,  questions  exist  relative  to 
the  achievable  height  accuracy  and  total  system  cost  and  complexity.  These  and  other  problem  areas 
are  discussed  later. 


PERFORMANCE  EVALUATION 

Each  of  the  above  systems  has  features  that  are  attractive  either  from  an  implementation  or  a 
system  perspective.  However,  these  qualities  are  meaningless  unless  the  design  can  provide  the 
required  topographic  precision  with  the  off-nadir  beam.  This  section  develops  a  detailed  analysis  of 
the  potential  performance  of  each  of  these  systems.  The  baseline  parameters  for  this  evaluation  are 
Altitude:  h  =  800  km 
One-sided  beam  offset:  50  km 
Equivalent  pointing  angle:  £  =  3.58° 

Frequency:  13.5  or  35  GHz 
Antenna  aperture  efficiency:  56% 

A  56%  aperture  efficiency  is  somewhat  conservative,  but  it  ensures  that  the  results  are  realistic 
rather  than  optimistic.  Note  that  only  the  outermost  beam  position  is  considered,  since  all  of  the  sys¬ 
tems  in  this  study  achieve  their  lowest  precision  at  this  outer  beam  position. 

Pure  Multiple-Beam  System 


A  key  parameter  in  estimating  the  precision  of  this  system  is  the  shape  of  the  average  return 
waveform.  This  waveform  is  the  convolution  of  three  functions:  the  flat-surface  impulse  response, 
the  system  point  target  response,  and  the  probability  density  function  of  the  surface  roughness.  When 
the  flat-surface  impulse  response  is  offset  from  nadir  (as  in  the  case  of  offset  beams)  and  its  time 
spread  (width)  is  considerably  larger  than  that  of  the  other  two  functions,  the  convolution  reduces  to 
the  product  of  the  flat-surface  impulse  response  and  the  combined  width  of  the  other  two  functions. 
When  the  system  point  target  response  and  the  surface  roughness  are  both  Gaussian  shaped,  the  com¬ 
bined  width  is  the  root  sum  square  of  the  widths  of  each  of  the  individual  functions.  Combining  this 
result  with  the  asymptotic  form  for  the  flat-surface  impulse  response  yields  [8] 


PA  r)  = 


nP-rG^V^mcT,/!) 

(4x)3ZTr3 


2t 


a(  1  t-  2b la) 


1/2 


exp 


-4 


1 


7(1  +  (2) 


[sin  £  —  €  cos  £]*■ 


O) 


where 

Pr  is  received  power, 

t)  is  pulse  compression  ratio, 

PT  is  peak  transmitted  power, 

G0  is  peak  boresight  antenna  gain, 

X  is  radar  wavelength. 

a°(i)  is  surface  cross  section  per  unit  scattering  area  at  the  angle  </- 
and  tan  \p  =  (cT/h)'/2  , 
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Lp  is  propagation  loss, 

$B  is  antenna  half-power  beamwidth, 

£  is  angle  between  the  antenna  boresight  and  nadir, 

Tc  is  composite  width  (time  spread)  of  the  system  point  target  response  and 
the  PDF  of  the  surface  roughness, 
h  is  altitude, 

c  is  speed  of  light, 

e  ~  ( cr/h)u 2, 

r  =  t  -  2h/c  (i  =  0  corresponds  to  the  time  when  the  pulse  is 

tiansmitted  by  the  radar), 

7  =  2.9  sin2(0fl/2), 

4  e  sin2£ 

a  =  - t~, 

7  0  +  *2) 

,  4  e2  sin2£ 

b  =  - 7-,  and 

7  (1  +  f2) 

r  =  h(\  +  e2)1'2. 

This  waveform  is  very  nearly  Gaussian  shaped,  as  it  should  be,  because  the  antenna  pattern  is 
assumed  to  be  Gaussian.  The  peak  of  the  return  power  waveform  occurs  where 


and  for  the  parameters  given  previously, 


r* i  .  2  y 

r0  =  —  tan  £, 
c 


tq  ~  10.44  (is. 


This  is  the  time  of  occurrence  of  the  peak  in  the  mean  or  average  return  and  corresponds  to  the  beam 
pointing  50  km  to  one  side  of  nadir.  The  shape  of  the  return  is  critical  in  designing  the  type  of  range 
tracker  used  to  lock-on  to  and  follow  this  mean  return.  The  design  of  one  type  of  tracker  is  discussed 
next;  however,  before  dealing  with  the  tracker  design,  it  is  convenient  to  derive  the  time  expanse  of 
the  waveform  given  by  Eq.  (1). 

First,  e  and  t  are  centered  on  the  peak  of  the  return  waveform,  i.e., 

«  =  €0  +  5e,  t  =  t0  +  6t,  (2) 

and  €0  =  ( CT0/h)xa .  Because  e  =  (cr/h)w2  and  6t/t0  «  1, 


e0  +  = 


1  + 

2  Tq 


6e  =  \  6r  ~~~ 
2  Ht0 
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Substituting  these  results  into  the  exponential  part  of  Eq.  (1)  yields 


exp 


-4 


7(1  +  e2) 


[sin  £  -  £  cos  £]2 


~  exp 


4_  cos2£ 
7  4 


r 

c 


(St)' 


The  value  of  6r  that  causes  this  exponential  function  to  fall  to  one-half  its  value  at  6r  =  0  is 


-  ln(0.5) 

The  total  half-power  width  is  2 5r  or 

2  h 


s 

s- 

K* 

J 

sin2£  1 

U2  J 

cos6£  j 

=  (6t)2. 


2St  =  —  tan  £  [-7  In  (0.5)/cos4£]1/2. 

c 


(5) 


(6) 


Substituting  from  the  relationship  between  7  and  the  ha’,  i- power  beamwidth  yields 

25 r  =  —  sin  (0fl/2)(1.42)  (7) 

c  cos2£ 

By  using  the  small  angle  aoproximations  for  both  £  and  »•©•» 

sin  (0Bi 2)  »  9b/ 2, 

tan  £  =  £, 

cos  £  ~  1, 

and  noting  that  for  a  56%  efficient  circular  aperture  [16], 

1.47X/D,  (8) 

where  D  is  the  diameter  of  the  aperture,  Eq.  (7)  becomes 

26r  =  2.08  (9) 

Df 

where  /is  the  frequency  of  the  radar.  The  corresponding  time  spread  Te  between  the  e-1  values  is 
Te  =  2.5  (/1  £//>/).  Figure  5  is  a  plot  the  half-power  width  for  the  present  case.  Equation  (9)  shows 
that  the  width  of  the  average  return  is  directly  proportional  to  the  altitude  and  pointing  angle  of  the 
radar’s  antenna  and  is  inversely  proportional  to  the  antenna  aperture  diameter  and  the  radar’s  fre¬ 
quency.  Having  both  the  average  return’s  shape  and  width,  an  expression  can  be  derived  for  the 
error  involved  in  making  a  waveform-based  measurement  of  the  off-nadir  range.  Figure  6  shows  a 
model  receiver  system  for  computing  the  rms  height  error. 
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Fig.  5  —  Plot  of  the  half-power  width  of  the  average  return  power  waveform 
for  a  56%  efficient  aperture  antenna 
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Fig.  6  —  Model  receive!  system  used  to  compute  rms  height  error 
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The  received  signal,  prior  to  detection,  will  comprise  a  zero-mean  narrow-band  Gaussian  pro¬ 
cess.  The  mean  square  value  (or  the  variance  <j2(/))  of  this  signal  is  equivalent  (within  a  scaling  fac¬ 
tor)  to  the  average  waveform  Pr(r)  discussed  previously.  If  the  additive  noise  variance  of  the 
receiver  is  denoted  by  a2,  the  signal  entering  the  square-law  detector  has  the  total  variance  equal  to 
a 2  +  a2  and  a  two-sided  (IF)  signal  bandwidth  equal  to  2/7^,  where  Tp  (~3  ns)  is  the  transmitted 
signal  compressed  pulse  length.  Let  Te  be  the  time  spread  of  Pr{r)  between  the  e~‘  points.  A 
theoretical  estimate  of  the  rms  height  error  for  a  split-gate  tracker  can  now  be  derived.  As  discussed 
in  Davenport  and  Root  [17],  for  a  Gaussian  input  signal  plus  noise  as  +  a„,  the  square-law-detected, 
zonal  filtered  signal  is  exponentially  distributed  with  the  variance 


=  a2  [jt; 


where  a  is  the  detection  scale  factor,  or 


oy  =  aa2 


1  + 


1 


SNR  SNR2 


1/2 


where  the  signal-to-noise  ratio  (SNR)  is  defined  at  the  peak  value  of  a2,  i.e.,  (o2snax/a2).  Replacing 
aa 2  by  the  range  spread  function  Pr(r)  and  considering  the  averaging  property  of  the  postdetection 
and  signal-averaging  filters,  we  obtain 


Pr(T) 


Oy  = 


Te 


1  + 


SNR 


SNR2 


1/2 


(10) 


where  (\/Tp)  and  (1/7,)  are  the  bandwidths  of  the  postdetection  filter  and  the  signal-averaging  filter, 
respectively.  Two  filters  are  discussed  for  conceptual  purposes  only;  in  an  actual  system  only  one 
composite  filter  would  be  needed.  Statistically  independent  samples  spaced  by  Tp  over  the  time 
expanse  Te  have  been  assumed.  The  factor  1/2  arises  from  replacing  the  idealized  zonal  filter  with  a 
realizable  filter  [18,19]. 


Next  we  use  the  waveform  model 


Pr(r)  =  e 


r 

T3 

2(T./2)\ 


and  assume  a  split-gate  tracker  that  tracks  the  Pr(r)  waveform  at  Pr(r)  =  e 
At/AP r  at  Pr(t)  =  e~112  is 


(ID 

Note  that  the  slope 


At 


Te 


A Pr  2e  ~0/2) 


We  next  determine  aT  =  ay  (At/ A Pr),  the  range  uncertainty  in  time  units,  as 


a'  =  7 


T  T 

tp  if 


>+^+  1 


SNR  SNR2 


1/2 
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To  convert  ar  into  distance  units,  we  use  r  =  (ct/2).  Furthermore,  considering  the  averaging  of  a 
number  of  statistically  independent  return  waveforms  over  a  time  period  TA ,  we  arrive  at  the  final 
result  for  range  uncertainty  ah  as 


Tp  Te 

(PRF)  (Ta) 


r 

i  +  -L- 

SNR 


+ 


1 

SNR2 


1/2 


(12) 


These  results  can  be  compared  with  the  Walsh  simulation  results  [20],  For  the  aircraft  case  with 
PRF  =  200,  Ta  =  0.05  s,  Tp  =  5  ns,  and  Te  =  33  ns,  Eq.  (12)  (for  SNR  =  100)  gives  ah  =  15.2 
cm.  The  simulation  yields  a  value  of  15  cm  [20], 

The  previous  discussion  has  tacitly  assumed  a  symmetrical  aperture,  with  equal  beamwidths  in 
both  planes.  The  analysis  is  readily  generalized  to  cover  unequal  along-track  and  cross-track 
beamwidths.*  For  the  general  case,  ah  can  be  shown  to  be  proportional  to  the  square  root  of  the 
cross-track  beamwidth  divided  by  the  along-track  beamwidth.  Thus  minimizing  ah  requires  minimiz¬ 
ing  the  cross-track  beamwidths  and  maximizing  the  along-track  beamwidths  (subject  to  resolution  con¬ 
straints). 


The  range  noise  expression  depends  on  a  number  of  parameters  that,  within  engineering  limita¬ 
tions,  can  be  arbitrarily  selected.  For  the  PRF,  however,  it  is  desirable  to  use  as  large  a  value  as 
possible.  The  largest  value  of  PRF  is  that  which  ensures  decorrelated  returns.  Larger  values  will 
reduce  the  effects  of  receiver  thermal  noise,  but  they  will  not  reduce  the  effects  of  the  inherent  ran¬ 
domness  of  the  return  power.  Since  thermal  noise  is  not  the  primary  problem  here,  it  does  no  good 
to  use  a  PRF  that  is  larger  than  the  one  for  decorrelated  returns. 

The  minimum  distance  necessary  for  the  radar  to  travel  for  decorrelation  may  be  accurately 
estimated  from  the  Van  Cittert-Zemike  theorem  [21J  as  modified  by  Walsh  in  Ref.  22  to  account  for 
the  two-way  radar  path.  In  particular,  this  distance  /  is  given  by 

l  =  0.3  h\/p,  (15) 

where  p  is  the  radius  of  the  illuminated  spot  on  the  surface.  For  the  case  at  hand  where  very  narrow 
beamwidths  are  considered. 


p  =  hdg/2. 


so  that 


(16) 


/  =  O.6X/0b. 


(17) 


Using  Eq.  (8)  to  relate  0B  and  the  antenna  aperture  diameter  yields 

t  =  0.408D.  (18) 


•Throughout  this  report,  we  assume  that  only  cross-track  squinted  beams  are  present;  other  geometries  may  be  treated  by  direct  analogy. 
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Since  the  PRF  is  given  by 


PRF  -  Vs/i, 


(19) 


therefore 


PRF  =  2.45  Vs/D, 


where  Vs  is  the  horizontal  velocity  of  the  altimeter.  Figure  7  is  a  plot  of  this  relation  for  the  present 
case.  Substituting  this  result  into  the  relationship  for  the  range  noise  yields 


0^  =  0. 126  c 


hijTc 

JsfTA 


1 

SNR2 


y  y 


1/2 


(20) 


ANTENNA  DIAMETER  (m) 

Fig.  7  —  Maximum  PRF  for  decorrelation  of  the  returns 
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This  result  predicts  higher  range  noise  with  altitude,  pointing  angle,  and  composite  pulse  width. 
Conversely,  a  lower  range  noise  can  be  expected  from  increasing  the  altimeter’s  horizontal  velocity, 
its  frequency,  the  averaging  time,  and  SNR.  Of  particular  interest  here  is  the  fact  that  the  range  noise 
is  independent  of  the  antenna  diameter  for  SNR  >  >  2.  This  is  a  bit  deceptive,  however,  since  it 
assumes  a  maximum  PRF.  That  is,  since  the  width  of  the  return  is  proportional  to  both  the  antenna 
diameter  and  the  maximum  PRF,  these  two  dependencies  cancel  because  they  appear  as  a  ratio  in  the 
expression  for  ah.  Assuming  that 

Tc  =  3.125  ns, 

Vs  =  7  km/s, 

Ta  =  Is, 

and  that  SNR  >  >  2,  Eq.  (20)  yields  the  following  range  noise  values  for  ah: 

f  4.9  cm  (13.5  GHz)  nu 

°h  ~  1  2.9  cm  (36  GHz)  ‘ 

If  constraints  dictate  that  the  PRF  cannot  be  the  maximum  value  for  a  given  antenna  diameter 
(as  shown  in  Fig.  7),  range  noise  will  increase.  Figure  8  shows  how  the  range  noise  increases  as  the 
PRF  decreases  (contracts)  from  its  maximum  value  for  a  frequency  of  13.5  GHz.  Figure  9  is  the 
same  plot  for  36  GHz.  These  curves  show  that  a  reduction  in  PRF  from  its  maximum  value,  by  a 
factor  a,  increases  the  range  noise  by  a1/2. 

Figures  10  and  11  are  plots  of  SNR  for  13.5  and  36  GHz,  respectively,  as  a  function  of  antenna 
diameter.  In  these  plots,  GEOSAT  altimeter  parameters  [23]  are  used: 

tj  =  32,768, 

PT  =  20  W, 
a0 (3. 58°)  =  10  dB, 

Receiver  NF  is  5  dB, 


Bandwidth  is  320  MHz. 


It  is  clear  from  these  plots  that  a  low  SNR  is  not  a  significant  problem. 
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Fig.  8  —  Increase  in  range  noise  with  decrease  or  contraction  of  the  PRF 
from  its  maximum  value  for  13.5  GHz 
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PRF  CONTRACTION 

Fig.  9  —  Increase  in  range  noise  with  decrease  or  contraction  of  the  PRF 
from  its  maximum  value  for  36  GHz 


16 


SNR  (dB) 


NRL  REPORT  9229 


We  next  consider  the  spatial  resolution  of  this  system.  For  mesoscale  opographic  mapping, 
some  degree  of  correspondence  should  be  provided  between  time  averaging  and  beamwidth  spatial 
averaging.  Spatial  averaging  may  be  incorporated  into  the  above  formulation  as  follows.  For  a 
Guassian  antenna  pattern  of  the  form  exp  (-B2 /lo2),  the  two-way  beamwidth  in  a-units  can  be  shown 
to  be 
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Using  Eq.  (8)  gives 


Ob 

or  for  13.5  GHz, 


D  <4.12  m.  (24) 

Figure  10  shows  that  the  SNR  for  this  antenna  diameter  is  —  30  dB,  so  a  low  SNR  is  not  a  problem. 
A  4-m  antenna  diameter  is  significant  from  another  prospective:  it  is  near  to  the  maximum  diameter 
antenna  that  a  Titan  HI  shroud  will  cover. 

Thus  if  a  4-m  antenna  diameter  is  used  along  with  the  maximum  PRF,  it  is  possible  to: 

•  achieve  a  <5-cm  range  noise  at  13.5  GHz, 

•  use  well-established  antenna  technology  to  generate  the  offset  beams,  and 

•  use  an  antenna  that  will  not  exceed  existing  launch  vehicle  capacity. 

From  these  results,  it  is  obvious  that  the  pure  multiple-beam  antenna  system  is  a  viable  means 
for  achieving  relatively  precise  off-nadir  altimetry. 

At  36  GHz,  the  propagation  considerations  are  that  ionospheric  effects  are  negligible;  however, 
precipitation  losses  are  significant.  A  precipitation  rate  of  12.5  mm/h  corresponds  to  an  attenuation 
of  =4  dB/km  [25].  For  a  storm  with  a  freezing  level  of  10  km,  the  two-way  attenuation  would  be  80 
dB.  The  corresponding  figure  for  a  rate  rf  2.5  mm/h  is  10  dB.  Thus  a  36  GHz  system  will  probably 
be  limited  to  rainfall  rates  of  <5  mm/h  fi-r  a  low-power  system  with  a  solid-state  transmitter.  A  final 
decision  on  operating  frequency  should  be  based  on  a  detailed  analysis  of  the  percent  coverage  loss 
that  this  precipitation  level  represents. 

Interferometric  Approach 

Figure  12  shows  the  interferometric  approach.  Two  antennas  are  separated  by  a  distance  d  in 
the  cross-track  direction  and  appropriately  phased  to  shift  their  composite  beam  off-nadir  by  an  angle 
£.  Since  the  desired  one-sided  cross-track  distance  is  50  km,  the  shift  angle  is  3.58°.  To  make  the 
footprint  on  the  surface  as  small  as  possible,  d  is  taken  to  be  many  wavelengths. 

The  derivation  of  the  flat-surface  impulse  response  for  this  configuration  is  somewhat  involved 
and  so  it  will  only  be  highlighted  here.  The  composite  antenna  pattern  is  a  product  of  the  pattern  of 
the  individual  antennas  and  the  so-called  array  factor.  Since  this  array  performs  best  if  its  boresight 
is  shifted  to  the  center  of  the  desired  spot  on  the  surface,  this  is  the  case  that  was  actually  analyzed. 
That  is,  slightly  better  performance  can  be  achieved  with  the  array  if  its  principal  beam  is  mechani¬ 
cally  rather  than  electronically  shifted  to  the  desired  spot  by  phasing  one  antenna  relative  to  the  other. 
Another  reason  for  doing  this  is  to  develop  a  general  expression  to  deal  with  array  pointing  errors. 

The  resulting  expression  for  the  composite  antenna  pattern  is  a  rather  complicated  function  of 
geometry.  To  simplify  the  azimuthal  integration  in  the  expression  for  the  flat-surface  impulse 
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Fig.  12  —  Interferometric  approach  to  forming  a  narrow  off-nadir  beam. 
The  interferometer  baseline  is  d. 


» 
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response,  certain  large  parameters  are  used  in  the  system  description.  However,  unlike  the  asymp¬ 
totic  evaluation  presented  in  Ref.  8,  the  proper  technique  to  evaluate  tne  integral  must  be  selected. 
That  is,  if  the  pattern  of  each  antenna  element  is  dominant,  Laplace’s  method  should  be  used.  How¬ 
ever,  if  the  array  factor  is  dominant,  stationary  phase  should  be  the  method  of  approximation.  A 
decision  can  be  based  on  the  relative  magnitude  of  k0d  compared  to  (4/y)  sin  2£.  The  former 
represents  the  dominance  of  the  array  factor  while  the  latter  represents  the  individual  antenna  element. 
Thus,  if 

kQd  >  —  sin  2i  »  1,  (25) 

7 

where  ko  is  the  radar  wavenumber  (k0  =  2t/X),  then  the  integration  can  be  accomplished  by  using 
the  stationary  phase  approximation.  Assuming  that  the  antenna  elements  comprise  56%  efficient  para¬ 
bolic  reflectors,  this  expression  can  be  reduced  to 

d  >  0.81  (£Z>2/X),  (26) 
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where  D  is  the  diameter  of  the  antenna  element.  If  the  equality  is  reversed,  Laplace’s  method  should 
be  used.  As  noted  in  Eq.  (25),  either  or  both  of  these  parameters  must  be  large  compared  to  unity, 
because  only  then  can  the  asymptotic  integration  techniques  be  considered. 


For  the  situation  described  in  Ref.  12,  the  inequality  described  in  Eq.  (26)  was  indeed  satisfied 
so  that  the  approximate  integration  was  accomplished  with  the  aid  of  stationary  phase.  This  led  to  the 
following  approximate  expression  for  the  flat-surface  impulse  response  when  using  the  two-element 
interferometer: 


P fs(t) 
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(27a) 


(27b) 


This  result  can  be  converted  to  the  actual  average  return  waveform  by  simply  multiplying  Eq.  (27)  by 
the  composite  pulse  width  Tc ,  provided  that  the  lobes  in  Eq.  (27)  are  mucii  wider  than  the  composite 
pulse  width.  The  angle  \p  is  given  by 


tan  \p  =  ( cr/h )1/2.  (28) 

Reference  12  provides  a  set  of  baseline  parameters  for  a  possible  spacebome  application;  these 
parameter  values  are 

d  =  10  m,  D  =  1.5m,  and  X  =  2.22  cm  (13.5  GHz). 
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Substituting  these  parameters  into  Eq.  (26)  yields  5.2  for  the  right-hand  side  and  lO  for  the  left-hand 
side;  thus,  the  inequality  is  indeed  satisfied  and  Eq.  (27)  can  be  used.  Figure  13  is  a  plot  of  the  flat- 
surface  impulse  response  for  the  10-m  interferometer.  The  lobes  are  due  primarily  to  the  second 
cosine  term  in  Eq.  (27).  The  peak  in  the  impulse  response  does  not  occur  exactly  at 

tan  £  =  (cT0/h)1'2 

because  of  the  x/4  in  the  argument  of  the  cosine  terms  in  Eq.  (27).  The  slight  although  noticeable 
asymmetry  in  the  shapes  of  the  lobes  is  due  to  the  two  cosine  terms  in  Eq.  (27)  Having  different 
periods.  The  failure  of  the  impulse  response  to  go  to  zero  between  the  peaks  in  the  lobe  structure  is 
due  to  the  nonzero  beamwidth  of  the  interferometer  elements  in  the  along-track  direction.  The  slight 
decrease  in  the  amplitudes  of  the  lobes  on  either  side  of  the  maximum  amplitude  lobe  is  due  to  the 
effect  of  the  1.5-m  interferometer  element. 


H  —  800  km  {-3.57*  ANTENNA  DIA  -  1.5  m  '.-I0.42w« 

f  —  50  km 


RELATIVE  DELAY  TIME  (ns) 

Fig.  13  —  Flat-surface  impulse  response  for  the  10-m  interferometer 


Figure  13  clearly  shows  the  major  drawback  to  the  interferometer  approach:  too  many  lobes  are 
very  similar  in  appearance.  From  a  measurement  standpoint,  it  would  be  impossible  to  tell  a  priori 
which  lobe  represents  the  return  from  the  desired  cel!  position  on  the  earth.  Designing  a  tracker  to 
“look  inside”  a  given  time  window  would  be  very  risky  because  of  the  possibility  that  the  tracker 
would  jump  to  the  adjacent  late  or  early  lobes.  This  basic  situation  is  further  complicated  by 
pointing-angle  uncertainty.  That  is,  if  errors  in  pointing  angle  and  direction  are  included,  there  may 
not  even  be  a  lobe  which  is  maximum  near  the  desired  spot  on  the  earth.  For  example,  it  may  be  that 
a  local  null  appears  at  r0  rather  than  a  local  maximum  (see  Fig.  13).  There  is  no  way  to  compensate 
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for  this  effect  unless  the  pointing  angle  and  direction  are  known  a  priori.  Since  the  pointing  direction 
is  very  difficult  to  measure,  this  is  a  major  drawback  to  the  interferometric  approach. 

A  possible  solution  to  this  lobe  identification  problem  is  to  increase  the  size  of  the  array-element 
antenna.  This  will  leave  the  main  lobe  unchanged  but  will  attenuate  the  adjacent  lobes  and,  thus, 
make  it  easier  to  track  the  main  lobe.  Increasing  the  size  of  the  array  element  means  increasing  D  in 
Eq.  (26);  this  will  lead  to  a  reversal  of  the  inequality  in  Eq.  (26)  and,  hence,  invalidate  Eq.  (27). 
Rather  than  redo  the  derivation  of  Eq.  (27)  when  the  inequality  in  Eq.  (26)  is  reversed,  there  is  an 
easier  method  of  estimating  the  effects  of  increasing  D.  By  increasing  D,  the  'obes  next  to  the  main 
lobe  near  tq  are  decreased.  The  total  time  between  the  lobe  just  before  the  peak  lobe  and  the  lobe 
just  after  it  is  —  1480  ns  (Fig.  13).  If  these  lobes  are  to  be  half  the  peak  lobe  amplitude,  we  can  use 
the  results  in  Fig.  5  to  estimate  the  diameter  of  the  antenna  element  that  will  accomplish  this.  A  time 
spread  of  1380  ns  requires  an  antenna  of  5.2-m  diameter  at  13.5  GHz.  However,  this  is  larger  than 
the  single  antenna  required  by  the  pure  multiple-beam  system.  Thus,  no  real  advantage  to  the  inter¬ 
ferometer  system  exists  relative  to  the  multiple-beam  system. 

Finally,  there  is  yet  another  problem  with  the  interferometer  concept.  In  the  actual  system,  the 
array  will  not  be  physically  tilted  to  point  to  a  cross  track  spot  50  km  away  from  nadir.  Instead,  one 
array  element  will  be  phased  relative  to  the  other  to  electronically  shift  the  beam  to  the  desired  spot. 
This  phasing  process  is  frequency  dependent.  That  is,  at  another  frequency  the  array  will  point  in  a 
different  direction.  Use  of  a  3.125-ns  compressed  pulse  width  requires  a  minimum  bandwidth  of  320 
MHz  that  at  13.5  GHz  center  frequency  is  sufficiently  large  to  cause  a  spot  size  shift.  That  is,  the 
location  of  the  beam  spot  on  the  surface  will  change  over  the  bandwidth  of  the  transmitted  pulse.* 
This  will  have  the  net  effect  of  making  it  even  more  difficult  to  determine  exactly  where  on  the  sur¬ 
face  the  beam  is  located. 

In  summary,  the  interferometric  approach  has  a  major  shortcoming  in  the  number  of  very  simi¬ 
lar  lobes  it  produces.  Its  sensitivity  to  pointing  errors,  signal  bandwidth,  and  a  10-m  phase  delay  path 
make  it  unsuitable  for  use  in  a  spacebome  application. 


Pseudo  SAR 

Synthetic  aperture  radar  (SAR)  is  a  mature  technique  for  generating  radar  reflectivity  images; 
traditionally  it  uses  short  pulses  to  provide  high  cross-track  resolution  (in  conjunction  with  a  flat  or 
spherical  surface  assumption)  and  Doppler  processing  to  provide  a  narrow  along-track  synthesized 
beam.  The  slant  range  (or  phase)  of  a  sequence  of  returns  from  cross-track  distance  x  is 

r(y)  =  [( h  -  Ah)2  +  *2]1/2  »  h  -  Ah  +  (x2 /2h). 

For  unfocused  SAR,  the  maximum  synthetic  aperture  length  is  (h\),/2.  For  typical  satellite  alti¬ 
tudes  at  Ku-band  the  angular  resolution  is  on  the  order  of  100  m.  This  so-called  unfocused  SAR 
resolution  may  be  obtained  by  excluding  Doppler-shifted  returns  outside  the  requisite  spectral  range— 
the  Doppler  beam-sharpening  mode.  In  the  focused  mode,  higher  angular  resolution  is  achieved 


•This,  in  effect,  causes  the  nulls  to  partially  fill  in. 
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through  the  process  of  phase  correcting  the  sequential  returns  prior  to  summation.  In  either  case, 
along  rack  resolution  is  enhanced  by  using  SAR  techniques  in  conjunction  with  radar  altimetry 
methods.  Cross-track  resolution  is  more  difficult  to  achieve.  In  this  case,  use  of  short  pulse  length 
alone  cannot  be  used  to  improve  the  cross-track  resolution.  For  altimetry,  a  surface  model  (e.g., 
spherical  earth)  cannot  be  assumed  since  both  slant  range  and  angle-of-arrival  information  are  needed 
to  unambiguously  determine  local  height.  Use  of  conventional  SAR  in  this  context  amounts  to  using 
one  measurement  (ranging)  to  determine  two  unknowns  (slant  range  and  angle-of-arrival).  This  prob¬ 
lem  has  been  explored  by  NRL  and  others  in  recent  years  [14,15].  These  works  show  that  an  inter¬ 
ferometric  form  of  SAR  is  needed  to  provide  both  SAR  resolution  and  altimetric  data.  The  ambiguity 
problems,  previously  discussed  in  conjunction  with  interferometric  systems  also  apply  here.  Inter¬ 
ferometric  SAR  requires  that  two  SAR  systems  be  deployed  across  a  reasonably  long  baseline.  Two 
such  phase-coherent  SAR  systems  are  not  within  well-developed  technology;  they  may  be  considered 
to  be  future-generation  sensors  relative  to  the  category  under  discussion.  Such  techniques  would  pro¬ 
vide  contiguous,  wide-swath  coverage,  with  higher  development  costs.  Within  the  context  of  existing, 
proven  technology,  we  return  to  the  topic  of  altimetric/ SAR  processing  to  enhance  along-track  resolu¬ 
tion  only. 

Reference  14  provides  an  excellent  summary  of  altimetric/SAR  capability  in  the  along-track-only 
case.  As  discussed  therein,  resolution  in  the  orthogonal  coordinate  must  be  provided  by  a  real  aper¬ 
ture.  In  the  numerator  of  Eq.  (12),  Tp  represents  pulse  width;  in  the  denominator,  PRF  represents 
the  pulse  repetition  frequency.  Since  the  azimuthal  resolution  is  restricted  by  the  need  for  pulse-to- 
pulse  decorrelation,  the  azimuthal  (or  cross-track)  resolution  is  the  only  variable  parameter  in  this 
case.  Because  this  is  a  denominator  term,  for  the  ocean  scatter  with  a  homogeneous  statistics  case  it 
should  be  as  large  as  practicable  within  previously  discussed  resolution  constraints.  For  sea  ice,  land, 
or  continental  ice  sheet  mapping,  the  functional  behavior  is  far  more  complex.  In  such  cases,  the  sur¬ 
face  topography  must  be  characterized.  For  example,  with  such  surfaces,  large-body  scattering  is 
often  dominant  and  statistical  homogeneity  is  grossly  violated. 

Pointing  Angle  Considerations 

For  pulse-length-limited,  nadir-pointing  altimetry,  the  attitude  angle  can  be  adequately  estimated 
by  examining  the  trailing  edge  of  the  return  waveform.  For  beamwidth-limited,  multibeam  altimetry, 
the  pointing  problem  is  much  more  severe.  A  5-cm  height  accuracy  requires  microradian-level  atti¬ 
tude  knowledge.  References  10,  12,  and  20  discuss  a  number  of  approaches  to  the  problem,  such  as 
orbit-crossing  analysis  with  data  from  geoidally  quiet  areas.  A  technique  suggested  by  these  studies  is 
the  use  of  a  conventional  GEOSAT-type  altimeter  to  perform  nadir  ranging  for  comparison  with  a 
multibeam  nadir  measurement.  Either  along-track  and  cross-track  Doppler  measurements  or  a 
nonradar  sensor  wouid  be  needed  to  determine  the  Cartesian  components.  In  any  case,  the  pointing 
angle  is  considered  to  be  an  unresolved  problem  area. 

CONCLUSIONS 

The  primary  conclusion  of  this  report  is  that  a  pure  multiple-beam  system  has  the  potential  for 
providing  three  beams  with  50-km  separation  and  measurements  of  range  to  a  precision  of  better  than 
5  cm  from  an  altitude  of  800  km.  In  addition,  it  appears  that  these  measurements  can  be  achieved  by 
using  existing  antenna  technology.  Future  efforts  should  be  directed  toward  experimental  verification 
of  pointing-angle  determination  and  antenna  development  to  provide  the  desired  number  of  cross-track 
beams  without  mutual  interference.  Because  of  the  near-nadir  operation  of  this  system,  polarization 
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appears  to  be  one  way  of  achieving  adjacent  beam  isolation,  as  originally  discussed  in  Ref.  21.  A 
spherical  reflector  with  compensating  feeds,  such  as  the  Arecibo  concept,  should  also  be  investigated. 

The  interferometric  approach  does  not  appear  to  be  as  well  suited  to  space  application  as  origi¬ 
nally  envisioned.  The  multiple-lobing  structure  of  the  interferometer  makes  it  very  difficult  to  isolate 
and  track  any  one  particular  lobe.  If  attitude  errors  are  not  properly  accounted  for,  it  becomes  almost 
impossible  to  track  a  specific  beam.  If  the  size  of  the  antenna  elements  comprising  the  array  is 
increased  to  attenuate  the  unwanted  lobes,  these  antennas  must  be  made  larger  than  the  single  antenna 
in  the  pure  multiple-beam  approach.  In  addition,  the  interferometer  is  a  phased  array  that  is  sensitive 
to  the  frequency  spectrum  of  the  transmitted  pulse.  Furthermore,  the  spatial  resolution  is  more  than  a 
factor  of  2  poorer  than  the  pure  multiple-beam  system  considered  in  this  report.  Considering  all  of 
these  and  other  problems  leads  to  the  conclusion  that  the  interferometer  or  SAR  approach  is  not  viable 
for  spacebome  application. 

SUPPORTING  AIRCRAFT  EXPERIMENTATION 

A  number  of  questions  regarding  the  viability  of  the  multibeam  altimetry  can  be  resolved  in  the 
laboratory  and  on  airborne  platforms.  An  airborne  system  is  essential  for  obtaining  the  much-needed 
data  from  the  open  ocean  to  validate  the  theoretical  computations  and  predictions  and  to  develop  and 
evaluate  real-time  measurements  and  off-line  data-processing  algorithms. 

Two  multibeam  altimeter  systems  for  airborne  platforms  are  being  developed,  a  13.5  GHz  sys¬ 
tem  at  the  Naval  Research  Laboratory  and  a  36  GHz  system  at  NASA/Wallops.  Currently  a 
Memorandum  of  Understanding  exists  between  the  two  laboratories  to  share  resources  and  perform 
experiments  that  will  complement  each  other’s  research  effort. 

The  multibeam  system  under  development  at  NRL  has  the  following  characteristics: 


Number  of  beams  2  to  5 

Frequency  13.5  GHz 

Transmitted  peak  power  1  kW 

Bandwidth  700  MHz 

Pulse  width  3  to  10  ns 

PRF  0.1  to  10  kHz 

Maximum  off-nadir  pointing  12.5° 

Beamwidth  <2° 

Nominal  operating  altitude  3000  m 


The  planned  areas  of  experimentation  at  NRL  are: 

•  Verifiy  basic  scattering  physics  relative  to  mean  value  waveform  shapes  and  fluctuation 
statistics 

•  Measure  the  sensitivity  of  waveform  to  sea  state  changes  as  pointing  angle  is  varied  from  0° 
to  12.5° 

•  Acquire  off-nadir  waveform  data  by  using  a  nadir  beam  tracker  to  control  the  off-nadir 
waveform  samplers;  use  these  off-nadir  data  to  evaluate  various  range  tracking  schemes 
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•  Determine  adjacent  beam  coupling 

•  Perform  pointing  angle  determination  studies 

•  Study  off-nadir  beam  gain  normalization  problems 

•  Observe  off-nadir  wind  speed  (o°) 

•  Gain  actual  experience  with  the  interferometer  ambiguity  problem 
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